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Abstract 

We discuss the ratio of the angular diameter distances from the source to the lens, D^s, 
and to the observer at present, D s , for various dark energy models. It is well known that the 
difference of D s s between the models is apparent and this quantity is used for the analysis of 
Type la supernovae. However we investigate the difference between the ratio of the angular 
diameter distances for a cosmological constant, (D<i s /D S ) A and that for other dark energy models, 
(Dds/Ds) other m this paper. It has been known that there is lens model degeneracy in using 
strong gravitational lensing. Thus, we investigate the model independent observable quantity, 
Einstein radius (#e), which is proportional to both D^s/Ds and velocity dispersion squared, a 2 . 
Dds/Ds values depend on the parameters of each dark energy model individually. However, 
(Dds/D s )^ — (D<i s /D s ) other for the various dark energy models, is well within the error of a v 
for most of the parameter spaces of the dark energy models. Thus, a single strong gravitational 
lensing by use of the Einstein radius may not be a proper method to investigate the property of 
dark energy. However, better understanding to the mass profile of clusters in the future or other 
methods related to arc statistics rather than the distances may be used for constraints on dark 
energy. 



1 Introduction 



Recent observations of high redshift Type la supernovae (SNe la) suggested that the expansion of 
the Universe is currently accelerating [1]. The cosmic microwave background (CMB) anisotropy 
data, indicating a spatially flat universe [2] containing a low value for the cold dark matter (CDM) 
density parameter [3], has confirmed that the Universe is dominantly made up of a component 
with negative pressure (dark energy) to make up the critical density today 

The cosmological constant and/or a quintessence field are the most commonly accepted can- 
didates for dark energy. Although the cosmological constant is simple and favored by current 
cosmological observations, there is 50 to 120 orders of magnitude discrepancy between theory and 
the measured value [4]. The quintessence, which might alleviate this problem is a dynamical scalar 
field leading to a time dependent equation of state parameter [5]. Various scalar field potentials 
for the quintessence have been investigated [6]. 

It is important to use various ways of checking for the existence of dark energy in addition 
to SNe la and CMB anisotropy constraints on dark energy. A number of other tests have been 
considered including "geometric" tests using standard cosmological methods (the galaxy cluster 
gas mass fraction [7], the location of CMB peaks [8, 9], the redshift-angular size [10], the strong 
gravitational lensing [11] - [14], fluctuations of the luminosity distance [15], etc.). 

The statistics of gravitational lensing of quasars (QSOs) by intervening galaxies can constrain 
on the cosmological constant [12]. Lensed images of distant galaxies in cluster, arcs or rings, may 
provide a bound on the equation of state parameter of dark energy [13]. While SNela is used 
to determine the luminosity distance itself, a gravitational lensing system can be used measure 
the ratio of angular diameter distances. Thus, the gravitational lensing system is regarded as an 
independent tool that complements SNe la as a probe of dark energy [14]. 

However, the lensing observations primarily depend on the parameters of lens models with 
minor dependence on cosmological parameters [16]. There is the lens model degeneracy in both 
the projected mass density profile and the circular velocity profile. It is shown that we need to 
measure the Einstein radius and the velocity dispersion within 0(1)% accuracy in order to put a 
constraint on ojde- 

In the gravitational lensing, one of the observable quantities without having any model depen- 
dence is the Einstein radius (9e), which is proportional to the velocity dispersion squared (cr^) and 
the ratio of the angular distances Dd s /D s , where Dd s is the distance from the lens to the source 
and D s is that from the source to the observer. With different values of cosmological parameters, 
we can have different values of D^/Ds, i.e. different values of Be- Thus, it might be used for 
probing the property of dark energy, ojde- However, there is an ambiguity in measuring a v . If the 
error of a v measurement is not within the differences of D^/Ds between different cosmological 
models, then we cannot distinguish the differences between models by measuring Be- 

This paper is organized as follows. In the next section we review the gravitational lensing 
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Figure 1: A general lensing system. L is the center of the lens, and the line through L and the 
observer O is the optical axis. /3 is an unperturbed angular position of the source relative to that. 
a is the deflection angle of a light ray, thus an image of the source is observed at position 0. 
However all angles are very small, we can replace the real light ray (SI'O) by its approximation 
SIO. 

system with the basic equations used in lensing observations. We also briefly mention the most 
popular lens models. We review the various aspects of errors in modeling lens in Sec. 3. In Sec. 4, 
we check both the differences of D s and the differences of D^s/Ds between the cosmological 
constant and other dark energy models. Our conclusion is in the last section. 



2 Gravitational Lensing and Isothermal Galaxy Models 

Figure 1 shows a simple lensing system [17]. Consider the source sphere S s , i.e. a sphere with 
radius D s , centered at the observer O and the deflector sphere Sd with radius Dd, i.e. the distance 
to the center of the lens L. In addition, consider the observer sphere S Q where the source would 
have angular position (3 (bold faced characters represent 2-dimensional vectors in this manuscript) 
if the light rays from the source S were not influenced by the gravitational field of the deflector. 
However, since light rays are bent by the lens, the straight line SO is no longer a physical ray 
path. Rather, there are light rays which connect the source and the observer but they are curved 
near Sd- One such ray SI'O is drawn, together with its approximation SIO, consisting of the two 
asymptotes of the real ray. The separation of the light ray from the optical axis, LI, is defined 
as the impact vector £ in the lens plane. The angle a between the two asymptotes SI and IO is 
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the deflection angle caused by the matter distribution L, 



& = 4Gj R2 X(Z')-±-±-d 2 Z' (2.1) 



In? 

where the integral is over the lens plane and £(£') is the surface matter density at position £' 
resulting from the projection of the volume mass distribution of the deflector onto the lens plane. 
This is valid when the gravitational field is weak, hence the deflection angle is small. The observer 
will thus see the source at the position on his sphere S . 

From the geometry of Fig.l, we can easily derive a relation between the source position de- 
scribed by the unlensed position angle (3 and the position of the images 6 = £/Dd of the source 

= O-^±&(t). (2.2) 

We can reexpress this equation by using the distance 77 = D s (5 from the source to the optical axis 
as 

V = ^£ - D dM® ■ ( 2 - 3 ) 

It is useful to rewrite lens equations (2.2) and (2.3) in dimensionless form by introducing a 
length scale £o i n the lens plane, which is called as the Einstein radius in the lens plane and a 
corresponding length scale 770 = £oAs/-Dd in the source plane. The Einstein radius in the lens 
plane is given by 

^ = ^AGM d .^^, (2.4) 

where Md is the mass of the lensing object. By use of the definition of the dimensionless vectors 
x = £/£o and y = 77/170 as well as the dimensionless surface mass density (convergence) k(x) = 
Ti(^qx) /Ti cr where the critical surface matter density S cr = D^D d ' ^ ne ^ ens equations are 
rewritten as 

/ x mix) , . 

y = x - a(x) = x — , (2.5) 

x 

where a is the scaled deflection angle, 

= ^atfbaO = U Sx'^x')^^ , (2.6) 

Dg^o 7T JR2 \x - x Y 

and m(x) is defined as 



m 



(x) = 2 f K(x')x'dx' . (2.7) 
Jo 

If the lensing object is a point mass, then the cross section a for strong lensing events is given by 
[18] 

a = < 2 • (2-8) 
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We need to specify the mass distribution S to solve the lens equation (2.5). One simple and 
analytic solution of a differential equation for the radial mass distribution is the singular isothermal 
sphere (SIS) [18]. Its surface mass density, projected onto the lens plane, is given by 

where a v is the velocity dispersion along the line of sight and the length scale £o = 47rcr 2 Dd j^ s d3 ■ 
With these we get the lens equation 

V = x-£- (2.10) 

|x| 

with the convergence 

<x) = ^-. (2.11) 

Due to its simplicity and the consistency with the matter distribution of galaxies, the SIS is 
frequently used as the gravitational lens model. In addition to SIS, a nonsingular isothermal 
sphere (NIS), the singular isothermal ellipsoid (SIE) and Navarro- Frenk- White (NFW) [19] are 
the most commonly used gravitational lens models. All of them are distinguished by their own 
surface mass densities as in the equation (2.9). 



3 Errors in Modeling Lens 

In this section we review systematic errors in the ratio Dd s /D s [14]. As we mentioned in the 
previous section, the various isothermal galaxy models are specified by their surface mass densities 
(£), which are the function of the velocity dispersion (a v ) along the line of sight. Isothermal 
ellipsoid models (SIE, NIE) have an additional necessary measurement, the ratio / of the minor 
axis to the major axis, £ = \J i\ + which is related to the ellipticity e by 

f = y/(l + e)/(l-e). (3.1) 
The lens equation gives an elliptical image of the Einstein ring with the minor and major axes 

0± = <Wl±e, (3.2) 

where 

9 E = ^a 2 v ^, (3.3) 

as given in Eq. (2.4). Thus D^/Dg can be determined by measuring a v , f, and 9e- 

Even though gravitational lensing has been used as a useful cosmological tool to probe the high 
redshift universe, there are several problems. Observational quantities depend on a lens model, 
which has inherent uncertainties in itself [20]. In most cases, we do not know the property of the 
lensing object in detail. Also, the light propagates through the local inhomogeneous spacetime 
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which deviates from the smoothed Robertson- Walker metric. Thus, even though the global pa- 
rameters such as the energy density parameters are fixed, the distance formula is not uniquely 
determined [21]. In spite of all these ambiguities, if dark energy dominates over the mass density, 
then the optical depth increases dramatically and its effect is much larger than the uncertainty 
arising from the problems in the formulation [22]. However in the following section we investigate 
the model independent observable quantity, Einstein radius. 



4 Angular Diameter Distances 

In this section, we investigate the ratio of the angular diameter distances from the lens to the 
sources (Dd s ) and those from the source to the present observer (D s ) for various dark energy 
models. We consider a spatially flat, homogeneous, and isotropic universe with radiation, matter, 
and dark energy. 

The angular diameter distance from the observer at present to the source, D s (0, z s ), is defined 

as 

d(0, z a ) 2 D(0,z a ) 

D s (0,z s ) = — = d L {l + z s ) =a s r s = 



Hq 1 { Zs dz 



(1+Z 8 



n£>(i + 2) 3(i + _g±5i_) + ng), e s/;(i-H*.)din(i + ^ 



(4.1) 



where di is the luminosity distance and z eq is the redshift when radiation and matter densities are 
equal. In the following numerical calculations, we use a eq = (1 + Zeq)^ 1 = 1/5510. The angular 
diameter distance depends on the present values of energy density contrast of matter (dm) and 
dark energy (O^e) and the equation of state of dark energy (cjde) as in Eqs. (4.1) and (4.2). In a 
flat universe, we have the constrain + f2j) E — 1 when we use the fact that fi£ 0) < tt$ , Q^. 
Especially, dark energy is the cosmological constant (A) when cj^e = — 1-0 and we call the model 
including cold dark matter and cosmological constant as ACDM model. 

However HqD s , itself may not be an interesting quantity in the study of gravitational lensing. 
Instead, we need to know the combination of angular diameter distances. We want to investigate 
the least model dependent observable quantity for our study and Einstein radius is one of the 
proper objects for this purpose. We can represent the angular diameter distance from the lens to 
the source as 

H^ 1 { Zs dz 



D ds (z d , z s ) = jt-j—t I i , (4.2) 

(1 + *.)■/* I n (0) {l ~ z)3 (! + ^) + n g), e 3/ (l^ DE )^ln(l +2 



where we assume z s > z&. 
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Figure 2: The ratio of angular diameter distances A const for the constant equation of state of 
dark energy (wde)- (a) Dependence on the present value of dark energy density contrast fi^E 
of A const for different present values of equation of state of dark energy Wpl, : —1.0 (solid line), 
—0.8 (dotted line), and —0.5 (dash-dotted line), respectively, (b) Dependence on the value of 

of A const for different values of : 0.85 (solid line), 0.73 (dotted line), and 0.60 (dash-dotted 
line), respectively. 



If we rewrite the Einstein radius, Eq. (3.3), then the ratio of the angular diameter distances 
D ds /D s is given by 

Dds 9e 



Rds 



D s Aval ' 

dz 



,/^ ) (l+ Z )3(l +T Ii±£ y )+^cxp(3 f z (l+^ DE )dln(l+ 2 ')) 

* " • ( 43 ) 



JZs 

^ ) (l+ Z )3(l +T Ii±fi y )+a^exp(3 £(l+u, DE )dln(l+ 2 ')) 

Thus, the error that contributes to this ratio due to the error from the velocity dispersion is 



Rds 



A ° V (4-4) 



We will not consider any detail of the lens model, which can arise additional errors. From various 
sources of the lens system, we can roughly say that the error of the velocity dispersion is about 
5% (Aa v /a v ~ 5%) [23]. That means we can distinguish the differences between models if the 
difference of the ratio of the angular diameter distances between models, ARds/Rds is at least 
5%. Now we will compare the differences of the ratios of angular diameter distances for different 
models. There was a work related to this approach [24]. However in the previous work, the 
equation of state parameter of dark energy was constrained as a constant. We consider more 
general cases with the stronger constraint on constant wdes. 
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Figure 3: A const for the constant cjde as a function of the location of the source (z s ). (a) De- 
pendence on the position of the source (z s ) of A const for different values of u>^, with the same 
notation for the label as the left panel of Fig. 2 when fij^ = 0.73. (b) Dependence on the position 
of lens (z,i) of A const for different present values of WpE, with the same notation as the left panel. 



4.1 R ds when u; DE = constant 

We assume that the stand cosmological parameters are ft^y^ = 0.73 and = —1.0 and 
we denote the ratio of the angular diameter distances D<is/D s with these parameters as R^ s . 
Figure 2 shows the difference of Rd s between the standard case and the various constant equation 
of state cases, A const = \R% S - R<£ nst \/R% s x 100(%). In the left panel of Fig. 2, we show the 
A const dependence on for different values of WpE, —1.0 (solid line), —0.8 (dotted line), and 
—0.5 (dash-dotted line). If there is an error in measuring fl^L within the range 0.51 < — 
0.87 (0.52 < < 0.89, 0.57 < ft£) for u§l = -1.0 (-0.8,-0.5), then A const is less than 5%. 
For example, if the measured cosmological parameters are = —0.5 and > 0.57, then 
we cannot distinguish this with the standard case, = -1.0 and = 0.73) within 5% by 

using Rd s , i.e. by measuring 9e- ^be a ^ so var i es f° r different values of wde, so = — 1.0 does 
not necessarily have smaller difference to the standard case as can be shown in the left panel of 
Fig. 2. When = —1.0, we can have the standard case if fij^ = 0.73. In this case A const = 0. 
This is shown in the left panel of Fig. 2 and it shows the consistency of figures. We show A const 
dependence on for different values of Q^m, 0.85 (solid line), 0.73 (dotted line), and 0.60 
(dash-dotted line) in the right panel of Fig. 2. Again if there is an error in measuring within 
the range u§l < -0.23 (-0.30, -0.44) for ng| = 0.85 (0.73,0.60), then A const is less than 5%. 
For example, if the measured cosmological parameters are = 0.85 and < —0.23, then we 
cannot distinguish this with the standard case within 5%. We can also see that A const = when 
fl^l = 0.73 and w£| = -1.0 in the right panel of Fig. 2. 

We show the dependence of A const on the position of source (z s ) in the left panel of Fig. 3 
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Figure 4: A p for the time varying equation of state of dark energy (wde = Wo + wi(l — a)), (a) 
Dependence on the value of ujq for different values of oj\ : 0.1 (solid line), 0.2 (dotted line), and 
0.3 (dash-dotted line), respectively, (b) Dependence on the value of lo\ of A p for different values 
of : —1.0 (solid line), —0.9 (dotted line), and —0.8 (dash-dotted line), respectively. 



for different values of —0.8 (solid line), —0.6 (dotted line), and —0.5 (dash-dotted line). We 
choose the position of lens as z& = 0.5 and = 0.73. If we choose z s = 3.0 for different values 
of wj^g, then we have the same values of A const of the dotted line of the right panel of Fig. 2. 
If the source position of lensing system is located z s > 0.64 (1.04) for w^e = —0.6 (—0.5), then 
A const is less than 5%. We cannot distinguish the difference of Rds s between the standard case and 
= —0.8 within 5% independent of the source position. In the right panel of Fig. 3, we show 
the dependence of A const on the position of lens (zd) when we choose z s = 3.0 and Qp E = 0.73 
with the same notation as the left panel of the figure for the different values of u>^- When we 
vary Zd from 0.1 to 1.5 for every values of (—0.8 < w^e — ~~ 0.5), A const s are well below 5%. 
To check the consistency between the figures, we can check A const values when we choose Zd = 0.5 
in the right panel of Fig. 2. This gives A const ~ 0.44 (1.5, 2.3) when cj£e = -0.8(-0.6, -0.5), 
which is equal to A const value when we choose z s = 3.0 in the left panel of the same figure. 

4.2 R ds when uj de = + <^i(l - a) 

Now we check the difference of the ratio of the angular diameter distances Dd s /D s between the 
standard case and one of the time varying lobe models, u>de = wo + u>i(l — a) [25]. We define 
A p = \Rd s — Rds\/R ds x 100(%) where R^ s means Rd s value when we use wde = <^o + <^i(l — a). 
We assume that there is no error in the measured value of ripg = 0.73 in this case. We show 
A p dependence on ujq value for different values of uj\, 0.3 (solid line), 0.2 (dotted line), and 0.1 
(dash-dotted line), respectively in the left panel of Fig. 4. A p is within 5% for uj\ = 0.3 (0.2,0.1) 
when ujq < —0.43 (—0.38, —0.34). For example, when uj\ = 0.3 and Uq < —0.43, then we cannot 
distinguish the cosmology with the cosmological constant from the cosmological model with time 
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Figure 5: A p . (a) Dependence on z s of A p for various sets of u)$ and lo\. (b) Dependence on Zd of 
A p for the same sets of loq and uj\ as in the left panel. 



varying dark energy wde = ojq + lo\{1 — a) within 5% by using Rd s . In the right panel of Fig. 4, we 

show the uj\ dependence of A p for different values of Uo, —1.0 (solid line), —0.9 (dotted line), and 

—0.8 (dash-dotted line). When u>\ varies from to 1 every A p value is within 5% for every value 

of ujq from —1.0 to —0.8. For the consistency check we can choose io\ = 0.3 for three different 

lines. We can check these values with the three different points of solid line of the left panel of 

the figure. When we choose loi = 0.3, A p = 0.450.76, and 1.2 for ujq = -1.0,-0.9, and -0.8 in 

both left and right panels of Fig. 4. 

We again check A p dependence on the position of z s and Zd for the various sets of ujq and io± 

values in figure 5. In the left panel of Fig. 5, we show A p dependence on z s . We choose Zd = 0.5 

and = 0.73 in this case. For the given sets of luq and u\ values, we cannot distinguish the 

cosmological model with the cosmological constant from the cosmological model with time varying 

dark energy lode = <^o + <^i(l — a) no matter where the source of lens is located from z s = 0.5 to 

6 within 5% by measuring Rds- We check Zd dependence on A p for the same sets of loq and tc\ 

in the right panel of Fig. 5. We choose z s = 3.0. Again none of case can be distinguished with 

cosmological constant case within 5%. 

4.3 R ds when u BE = -1 + Aa q c /3{a q + a q c ) 

We check another time varying dark energy model by using the parameter u>de = — l + 4aj?/3(a 9 + 
a q c ) in figure 6. We define A SL = \R^ S - R^ \/R^ 3 x 100(%) where R^ means R ds value when 
we use cjde = — 1 + 4a q /3(a q + a q ) [6]. The angular diameter distance of this parametrization is 
almost degenerate with that of the cosmological constant. Thus, the ratio difference of angular 
diameter distances between this parametrization and the cosmological constant is extremely small 
for the different choice of parameters. We choose fij^, = 0.73, z s = 3.0, and Zd = 0.5 in this 
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Figure 6: A for the varying equation of state of dark energy (wde = — 1 + Aa q c /?>{a q + a-?)), (a) 
Dependence on q of A SL for different values of x c '■ —2.00 (solid line), —2.64 (dotted line), and 
—3.00 (dash-dotted line), respectively, (b) Dependence on x c of A SL for different values of q : 1.0 
(solid line), 3.0 (dotted line), and 5.0 (dash-dotted line), respectively. 



figure. In the left panel of Fig. 6, we show the dependence of A on q for the different values of 
x c , —2.00 (solid line), —2.64 (dotted line), and —3.00 (dash-dotted line). For all of the cases, 

A SL 

are within 1% for different values of x c from q = 2 to 5. In the right panel of Fig. 6, we show A SL 
dependence on x c for different values of q, 1 (solid line), 3 (dotted line), and 5 (dash-dotted line). 
Except for q = 1 case, 

A SL 

are degenerate within 1%. Thus, we cannot tell any difference between 
the cosmological constant and the dark energy with the parameter cjde = — 1 + 4a^/3(a 9 + a q ) by 
measuring R^s within 5%. 

We again check A SL dependence on the position of z s and z^ for various sets of x c and q 
values in Fig. 7. In the left panel of Fig. 7, we show A SL dependence on z s . We choose Zd = 0.5 
and f2pg = 0.73 in this case. For given sets of (x c , q) values, (—2.00, q = 4) (solid line), (—2.64, 
q = 3) (dotted line), and (—3.00, q = 2) (dash-dotted line), we cannot distinguish the cosmological 
model with the cosmological constant from the cosmological model with time varying dark energy 
cl>de = — 1 + ^a q /2>{a q + a q c ) no matter where the source of lens is located from z s = 0.5 to 6 within 
1% by measuring Rd s - We check za dependence on A SL for the same sets of x c and q in the left 
panel of the figure. We choose z s = 3.0. Again none of the cases can be distinguished from the 
cosmological constant case within 1%. 



5 Conclusions 



We have investigated the angular diameter distances from the source to the present observer for 
the various dark energy models. The angular diameter distance is proportional to the luminosity 
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Figure 7: A SL . (a) Dependence on z s of A SL for various sets of x c and q values, (b) Dependence 
on Zd of A SL for the same sets of x c and q values as in the left panel. 



distance, which is used for the analysis of Type la supernovae. The difference of HqD s between 
models is apparent and measuring this quantity as to probe dark energy models might be proper 
for further investigation. In addition to SNela, the gravitational lensing has been known as the 
possible method to probe the property of dark energy. 

Most of the lensing analyses have the Einstein radius as a basic quantity. This is lens model 
independent observable quantity and might be a good probe for the property of dark energy. The 
Einstein radius is proportional to both the ratio of the angular diameter distances and the velocity 
dispersion squared. However, we have shown that there is degeneracy between different models 
for the value of the ratio of the angular diameter distances. If the error in measuring the velocity 
dispersion exceeds the difference of the ratio of the angular diameter distances of two different dark 
energy models, then we cannot distinguish the differences between different dark energy models 
by measuring the Einstein radius. In our analysis we have shown that this is the case for most of 
the parameter spaces of the dark energy models. Thus, a single strong gravitational lensing might 
not be a proper method to investigate the property of dark energy. However, better understanding 
to the mass profile of clusters in the future or other methods related to arc statistics rather than 
the distances may be still used for constraints on dark energy. 
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